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Specific heat C P (T) measurements of the heavy fermion 
normal-spinel structure compound LiV2C>4 were carried out 
using a heat-pulse calorimeter over the temperature T range 
from 1.2 to 108 K. The electronic specific heat C C (T) of 
LiV204 is extracted from the C P (T) data using the lattice 
contribution obtained for LiTi2 04, a superconductor with T c 
= 11.8 K. The electronic specific heat coefficient j(T) = 
C C {T)/T of LiV 2 4 is found to be 0.42 and 0.43J/molK 2 
at T = IK for two different high magnetic purity samples, 
respectively. f(T) decreases rapidly with increasing temper- 
ature from 4 to 30 K and then decreases much more slowly 
from 0.13J/molK 2 at 30 K to 0.08J/molK 2 at 108 K. The 
C e (T) of the first of the above two L1V2O4 samples is com- 
pared with theoretical predictions for the spin S = 1/2 Kondo 
model, a generic Fermi liquid model, and an antiferromagnet- 
ically coupled quantum-disordered metal. Each of these the- 
ories can adequately describe the T dependence of C e in the 
Fermi liquid regime at low (~ 1-10 K) temperatures, consis- 
tently yielding a large extrapolated 7(0) = 428(3) mj/mol K 2 . 
However, none of these theories describes C e (T) from ~ 10 K 
to 108 K. Our C C (T) data are also in severe disagreement 
with the magnetic specific heat of the spin S = 1/2 Heisen- 
berg model, calculated above ~ 40 K for the V sublattice 
of the spinel structure. Thermal expansion measurements of 
LiV204 were carried out from 4.4 to 297 K using a differen- 
tial capacitance dilatometer. Strong increases in the thermal 
expansion coefficient and Griineisen parameter F are found 
below ~ 20 K, confirming the results of Chmaissem et al. 
[Phys. Rev. Lett. 79, 4866 (1997)] obtained using neutron 
diffraction. We estimate T(0) ~ 11.4, which is intermediate 
between those of conventional metals and /-electron heavy 
fermion compounds. 



I. INTRODUCTION 

Heavy fermion (HF) and related intermediate valence 
(IV) behaviors are ubiquitous in metallic /-electron sys- 
tems containing lapthanide or actinide (= M) atoms with 
unstable valence!!! The HF materials are typically inter- 
metallic compounds containing Ce, Yb or U ions and are 
characterized at the lowest temperatures T by a large 
and nearly T-independent spin susceptibility \(T — > 
0) ~ 10~ 2 cm 3 /(mol M), and an extraordinary large 
nearly T-independent electronic specific heat coefficient 
7(T -f 0) ~ U/(mol M) K 2 , where 7 (T) = C e (T)/T 
and C e (T) is the electronic contribution to the mea- 



sured specific heat at constant pressure C P (T). Large 
quasiparticle effective masses m* of ~ 100-1000 elec- 
tron masses m c have been inferred from 7(0) for the 
HF compounds and smaller values for the IV materials. 
The normalised ratio of x(0) to 7(0), the Sommerfeld- 
Wilson ratio 2 ] Rw, is on the order of unity in HF and 
IV materials as in conventional metals, and is given by 
Rw = 7r2 fcix(0)/3Mcff7(0)' where fee is Boltzmann's con- 
stant and /Lt e ff is the effective magnetic moment of the 
Fermi liquid quasiparticlcs. For quasiparticles with (ef- 
fective) spin S — 1/2, one obtains 



Rw 



4^ 2 fcjx(0) 
3ff 2 ^7(0) 



(1) 



where g is the g-factor of the quasiparticles and hb is 
the Bohr magneton. Since Rw ~ 1 in many of the HF 
and IV compounds, x an( i Ce at low temperatures are 
both probing the same low-energy heavy quasiparticle 
spin excitations. With increasing T in the heaviest-mass 
systems, x{T) crosses over to local-moment behavior and 
7 decreases rapidly, on a temperature scale of ~ 0.3-30 K. 

Heavy fermion behaviors are not expected for d- 
electron compounds because of the much larger spa- 
tial extent of d orbitals than of / orbitals and the re- 
sulting stronger hybridization with conduction electron 
states. Recently, however, in collaboration with other re- 
searchers, we have documented HF behaviors, character- 
istic of those oithe heaviest mass /-electron HF systems, 
in the metallkp transition metal oxide compound LiV2C>4 
using C P (T)^ xiTAffa 7 Li and 51 V NMRjM muon 1 
relaxation (/kSR),ui3 and 4-295 K crystallography 
measurements- —Independent crystallogranityj-and x(T) 
measurementaJ'tJ and NMR measurementst3il3 were re- 
ported nearly simultaneously by other groups, with sim- 
ilar results. LiV204 has the face-centexed-cubic normal- 
spinel structure (space group Fd3m)]Ei and is formally 
a g? 1 - 5 system. The Li atoms occupy tetrahedral holes 
and the V atoms octahedral holes in a nearly cubic-close- 
packed oxygen sublattice, designated as Li[V2]C>4. The 
C C (T) is extraordinarily large for a transition metal com- 
pound, 7(1 K) w 0.42 J/molK?, decreasing rapidly with 
T to ~ 0.1 J/molK 2 at 30Kfl As discussed extensively 
in Refs. |and|, from - 50-100K to 400 K, x{T) shows 
a Curie- Weiss-like [x = C/(T — 6)] behavior correspond- 
ing to antiferromagnetically coupled (9 = —30 to — 60 K) 
vanadium local magnetic moments with S — 1/2 and 
g » 2, but static magnetic ordering does not occur above 
0.02 K in magnetically pure L1V2O4, and superconduc- 
tivity is not observed above 0.01 K. 
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To our knowledge, in addition to L1V2O4 the only 
other stoichiometric transition metal spinel-structure ox- 
ide which is metallic to low— temperatures is the normal- 
spinel compound LiTi204o1IZl In contrast to L1V2O4, 
this compound becomes superconducting at T c < 13.7 K 
(Refs. yQ and has a comparativel y j T- iiadependent 
and small x(T) from T c up to 300 kMI^] The re- 
sistivity of thin films at 15 K is (4.3-8.8) xl0~ 4 r2cm.B 
The spinel system Lii+^^-zCU with cation occuflpjez. 
LilLi^Tia-Jo^xists from x = to x = l/3f5BEZI 
for x — 1/3 f3 the oxidation state of the Ti is +4 
and the compound is a nonmagnetic insulator. A zero- 
temperature s up i er f jov iductor- insulator transition occurs 
at x ~ 0.1-0.2BBE3 

In this paper, we report the details of our C P (T) 
measurements on L1V2O4 and of the data analysis and 
theoretical modeling. We have now obtained data to 
108 K, which significantly- extends our previous high- 
temperature limit of 78K.Q We also present complemen- 
tary linear thermal expansion a(T) measurements on this 
compound from 4.4 to 297 K. We will assume that C P (T) 
can be separated into the sum of electronic and lattice 
contributions, 



the value at T = 0. Theoretical modeling of the C e (T) 
data for L1V2O4 is given in Sec. [v|. Since the electri- 
cal resistivity data for single crystals, of L1V2O4 indicate 
metallic behavior from 4K to 45OK,0we first discuss the 
Fermi liquid description of this compound and derive the 
effective mass and other paramete rs fo r the current car- 
riers at low temperatures in Sec. VA. This is followed 
by a more general discussion of the FL theory and its 
application to L1V2O4 at low T. In Sec. VB we compare 
the predictions of Zulicke and Milliso for a quantum- 
disordered antiferromagnetically coupled metal with our 
C e (T) results for LiV 2 4 . The isolated 5=1/2 impurity 
Kondo model predicts FL behavior at low temperatures 
and impurity local moment behavior at high tempera- 
tures. Precise predictions for the x(T) anc ^ Ce(T) have 
been made for this model, and we c omp are our C C (T) 
data with those predictions in Sec. |V C| . In Sec. VD 



we consider a local moment model in which the mag- 
netic specific heat of the B sublattice of the A[£> 2 ]C>4 
spinel structure for spins S — 1/2 and 5 = 1 per B ion 
is given by a high-temperature series expansion and the 
predictions compared with the C e (T) data for L1V2O4. 



C p (T) = C (T) + C lat (T) 



(2a) 



A summary and concluding remarks are given in Sec. VI 
Unless otherwise noted, a "mol" refers to a mole of for- 
mula units. 



C (T) = 7 (T)T . (2b) 

In Ref. |4|, we reported C P (T) measurements up to 108 K 
on Li4/ 3 Ti 5 / 3 04 which were used to estimate C lat (T) in 
L1V2O4 so that C C (T) could be extracted according to 
Eq. (pa|). In the present work, we report C P (T) up to 
108 K for LiTi204, compare these data with those for 
Li 4 / 3 Ti 5 / 3 04, and obtain therefrom what we believe to 
be a more reliable estimate of C lat (T)for LiV 2 4 . The 
experimental details are given in Sec. P. An overview of 
our C P (T) da ta for L1V2O4, LiTi2C>4 and Li 4 / 3 Ti 5 / 3 C)4 is 
given in Sec. Ill A . Detailed analyses of the data for the 
Lii+ 2; Ti2- : r04 compound s and comparisons with litera- 
ture data are given in Sec. Ill B| , in which we also estimate 
C lat (T) for L1V2O4. The C C (T) and electro nic entropy 
S C (T) for L1V2O4 are derived in Sec . |III C\ The a(T) 



measurements are presented in Sec. IV and compared 
with the C P (T) results and lattice parameter data versus 
temperature obtained from neutron diffraction measure- 
ments by Chmaissem et alB From the combined a(T) 
and C p (T) measurements on the same sample, we derive 
the Gruneisen parameter from 4.4 to 108 K and estimate 



II. EXPERIMENTAL DETAILS 

Polycrystalline L1V2O4 samples were prepared using 
conventional ceramic techniques described in detail else- 
where, where detailed sample characterizations and mag- 
netic susceptibility results and analyses are also given.cl A 
few of these results relevant to the present measurements, 
analyses and modeling are given in Table | 

Polycrystalline LiTi2C>4 and Li4/ 3 Ti 5 / 3 04 samples 
were synthesized using solid-state reaction techniques. L3 
TiC>2 (Johnson Matthey, 99.99 %) was dried under a pure 
oxygen stream at 900 °C before use. This was mixed with 
Li2C0 3 (Alfa, 99.999%) in an appropriate ratio to pro- 
duce either Li4/ 3 Ti 5 / 3 04 or a precursor "LiTiO^.s" for 
LiTi2C>4. The mixtures were then pressed into pellets 
and heated at 670 °C in an oxygen atomosphere for sa 1 
day. The weight loss due to release of CO2 was within 
0.04 wt.% of the theoretical value for LiTi02.5- How- 
ever, for Li 4 / 3 Ti 5 / 3 04 additional firings at higher tem- 
peratures (up to 800 °C), after being reground and repel- 
letized, were necessary. LiTi2C>4 was prepared by heating 



TABLE I. Lattice parameter ao and structural [/imp (Str)] and magnetic [/imp (Mag)] impurity concentrations for the LiV2C>4 
samples studied in this workE 



Sample No. 


2 


3 


4A 


5 


6 


Lattice parameter (A) 


8.23997(4) 


8.24100(15) 


8.24705(29) 


8.24347(25) 


8.23854(11) 


Impurity phase 


v 2 o 3 


pure 


v 2 o 3 


V 2 3 


v 3 o 5 


/imp (Str) (mol%) 


1.83 


< 1 


1.71 


< 1 


2.20 


/ m p (Mag) (mol 96) 


0.22(1) 


0.118(2) 


0.77(2) 


0.472(8) 


0.0113(6) 
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pressed pellets of a ground mixture of the LiTiC>2.5 pre- 
cursor and Ti2C>3 in an evacuated and sealed quartz tube 
at 700 °C for one week and then air-cooling. The Ti203 
was prepared by heating a mixture of TiC>2 and titanium 
metal powder (Johnson Matthey) at 1000 °C for one week 
in an evacuated and sealed quartz tube. 

Powder x-ray diffraction data were obtained using a 
Rigaku diffractometer (Cu Ka radiation) with a curved 
graphite crystal monochromator. Rietveld refinements 
of the data were carried— put using the program "Ri- 
etan 97 ('beta' version)" rB The x-ray data for our sam- 
ple of Li4/ 3 Ti 5 / 3 04 showed a nearly pure spinel phase 
with a trace of TiC>2 (rutile) impurity phase. The 
two-phase refinement, assuming the cation distribution 
Li[Li!/ 3 Ti5/3]04, yielded the lattice ao and oxygen u pa- 
rameters of the spinel phase 8.3589(3) A and 0.2625(3), 
respectively; the concentration of TiC>2 impurity phase 
was determined to be 1.3 mol%. The LiTi204 sample was 
nearly a single-phase spinel structure but with a trace 
of Ti 2 03 impurity. A two-phase Rietveld refinement 
assuming the normal-spinel cation distribution yielded 
the spinel phase parameters ao = 8.4033(4) A and u — 
0.2628(8), and the Ti 2 03 impurity phase concentration 
< 1 mol%. i-Qur crystal data at©, compared with those of 
Cava et alH and Dalton et alH in Table 0. 

The C P (T) measurements were done on samples from 
four different batches of LiV2 04 using a conventional 
heat-pulse calorimeter, with Apeizon-N grease providing 
contact between the sample and the copper tray.c3 Addi- 
tional C P (T) data were obtained up to 108 K on 0.88 g of 
the isostructural nonmagnetic insulator spinel compound 
Li4/ 3 Ti5/ 3 04, containing only maximally oxidized Ti +4 , 
and 3.09 g of the isostructural superconductor LiTi 2 04 
to obtain an estimate of the background lattice contribu- 
tion. A basic limitation on the accuracy of these C p data, 
except for LiV204 below 15 K, was the relatively small 
(and sample-dependent) ratios of the heat capacities of 
the samples to those associated with the tray (the 



"addenda"). For L1V2O4 sample 6, this ratio decreased 
from 40 near IK to 1.0 at 15 K to a relatively constant 
0.2 above 40 K. For the superconducting LiTi204 sam- 
ple, this ratio was 0.45 just above T c (= 11.8 K), and in- 
creased to 0.65 at 108 K. For the nonmagnetic insulator 
Li4/ 3 Ti 5 / 3 04 sample, this ratio varied from 0.03 to 0.12 
to 0.2 at 8, 20 and 108 K, respectively. These factors are 
important since small (±0.5%) systematic uncertainties 
in the addenda heat capacity can have differing effects 
on the C P (T) measured for the different samples, even 
though the precision of the raw heat capacity measure- 
ments (as determined from fits to the data) is better than 
0.25%. 

The linear thermal expansion coefficient of L1V2O4 
sample 6 was-,«ieasured using a differential capacitance 
dilatometer.ESEZI All data were taken isothermally (T 
constant to 0.001 K). The absolute accuracy of the mea- 
surements is estimated to be better than 1%. 



III. SPECIFIC HEAT MEASUREMENTS 
A. Overview 

An overview of our C P (T) measurements on L1V2O4 
sample 2, run 2 (1.26-78K), sample 6 (1.16-108 K), and 
LiTi2C>4 and Li4/ 3 Ti5/ 3 04 up to 108 K, is shown in plots 
of C* P (T) and C P (T)/T in Figs. 0(a) and (b), respec- 
tively. Our data for LiTi2 04 and Li4/ 3 Ti 5 /304 are gept 
erally in agreement with those of McCallum et aZ.cJ 
which cover the range up to ~ 25 K. For LiTi 2 04 above 
T c = 11.8 K (see below) and for Li 4 / 3 Ti5/ 3 04, one sees 
from Fig. |l|(a) a smooth monotonic increase in C p up to 
108 K. From Fig. 0(b), the C p of the nonmagnetic insula- 
tor Li4/ 3 Ti 5/ / 3 04 is smaller than that of metallic LiTi204 
up to ~ 25 K, is larger up to ~ 45 K and then becomes 
smaller again at higher temperatures. Since C c = in 
Li4/ 3 Ti 5 / 3 04 and C e (T) in LiTi204 cannot be negative, 



TABLE II. Characteristics of LiTi204 and 1^4/3^5/304 samples. Abbreviations: ao is the lattice parameter, u the oxygen 
parameter, 7 the electronic specific heat coefficient, 80 the zero-temperature Debye temperature, T c and AT C the supercon- 
ducting transition temperature and transition width, and AC P the specific heat jump at T c . 



ao 


11 


7 


0o 


T c 


AT C 


AC P /7T C 


Ref. 


(A) 




(mJ/molK 2 ) 


(K) 


(K) 


(K) 


(mJ/molK 2 ) 






LiTi2 04 


8.4033(4) 


0.2628(8) 


17.9(2) 


700(20) 


11.8 


^0.2 


1.75(3) 


This Work 


8.4033(1) 


0.26275(5) 














16 


8.41134(1) 


0.26260(4) 














17 


8.407 




21.4 


685 


11.7 


1.2 


1.59 




28 






22.0 


535 


12.4 


0.32 


1.57 




I 




0.26290(6) (300 K) 














21 




0.26261(5) (6K) 














21 


Li 4 /3Ti 5 / 3 04 


8.3589(3) 


0.2625(3) 





725(20) 








This Work 


8.35685(2) 


0.26263(3) 
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8.359 







610 










28 






0.05 
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FIG. 1. Overview of the molar specific heat C p (a) and 
Cp/T (b) vs. temperature T for LiV204 samples 2 (•) and 6 
(o) and the reference compounds LiTi2 04 (•), a metallic su- 
perconductor, and Li4/ 3 Ti 5/ /304 (o), a nonmagnetic insulator. 
The solid curves are polynomial fits to the data for L1V2O4 
sample 6 and LiTi2C>4. The dashed curve in (b) is the inferred 
normal state C p /T below T c for LiTi2 4 . 



it follows from Eq. (2a) that C lat (T) and hence the lattice 
dynamics are significantly different in LiTi204 compared 
with Li 4/ 3Ti 5/3 04. The data for LiV 2 4 in Fig. 0(b) 
are shifted upwards from the data for the Ti spinels, 
with a strong upturn in C P (T)/T below ~ 25 K. These 
data indicate a very large "f(T — > 0). Comparison of 
C p (T)/T for LiV 2 4 and LiTi 2 4 at the higher tem- 
peratures > 30 K indicate that a large j(T) persists in 
LiV 2 4 up to our maximum measurement temperature 
of 108 K. In the following, we begin our analyses with 
the data for the Lii +a; Ti 2 _ a ;0 4 compounds because we 
extract a lattice specific heat from these materials as a 
reference for LiV 2 4 . 



B. Lii+ x Ti2-i04 

In the present paper, our C P (T) data for LiTi 2 04 and 
Li4/ 3 Ti 5 / 3 04 are most important for determining the lat- 



tice contribution C lat (T) to C P {T) of LiV 2 4 . At low 
temperatures, the C P (T) of a conventional nonmagnetic, 
nonsuperconducting material iscB 



C P (T) = A\T + A 3 T 3 + A b T b + A 7 T 7 



(3) 



where 7 = A\ and [3 = A3. From Eqs. (|2j), the first term 
in Eq. (||) is C C (T), the second corresponds to the ideal 
Debye lattice contribution C lat (T — > 0) and the follow^ 
ing terms represent dispersion in the lattice properties— 
The zero-temperature Debye temperature 9q is given byE£ 
6> = (1.944 x 10 6 r/f3) 1 ^ 3 , where r is the number of 
atoms per formula unit (r = 7 here) and (3 is in units of 
mJ/molK 4 . Equation (^) suggests the commonly used 
plot of C p /T versus T 2 to obtain the parameters 7 and 
(3. Unfortunately, the very small heat capacity of the 
small Li4/ 3 Ti 5 / 3 04 sample and the occurrence of the su- 
perconducting transition in LiTi 2 04 at 11.8 K complicate 
the use of this relation to determine C lat (T) for these pre- 
sumably similar materials below ss 12 K. 

The C P (T)/T of LiTi 2 4 below 20 K is plotted ver- 
sus T and T 2 in Fig. ||(a) and (b), respectively. The 
superconducting transition at T c = 11.8 K is seen to be 
pronounced and very sharp (AT C < 0.2 K). The dotted 
line extrapolation of the normal state (T > 11.8 K) data 
to T = shown in Figs. |l](b) and [2] uses Eq. (||), equality 
of the superconducting and normal state entropy at T c , 
S^ll.SK) = 241(1) mJ/molK, and continuity considera- 
tions with C P (T)/T above T c , from which we also obtain 
estimates of 7 and f3. Although we cannot rule out a In- 
dependence of 7, we assume here that 7 is independent 
of T. While 7 [= 17.9(2) mJ/molK 2 ] appears to be quite 
insensitive to addenda uncertainties, 9 [= 700(20) K] 
is less well-defined. Our value for 7 is slightly smaller 
than the|-s4alues of 20-22 mJ/molK 2 reported earlier for 
LiTi 2 04p3€3 as shown in Table ||. From the measured 
superconducting state C P (T C ) = 684(2) mJ/molK and 
normal state C P (T C ) = 315(1) mJ/molK, the discontinu- 
ity in C p at T c is given by AC P /T C = 31.3(3) mJ/molK 2 , 
yielding AC p /-fT c = 1.75(3) which is slightly larger than 
previous estimates in Table |J. According to Eqs. (||), 
the lattice specific heat of LiTi 2 04 above T c is given by 
C lat (T) = C P (T) -7T. 

The C lat (T) derived for LiTi 2 4 below 12 K is consis- 
tent within experimental uncertainties with the measured 
C lat (T) of Li 4 / 3 Ti 5 / 3 04 in the same temperature range 
after accounting for the formula weight difference. The 
low-T CJT)/T = C ut {T)/T for Li 4/3 Ti 5/3 04 is plotted 
in Figs. U The 6 = 725(20) K found for Li 4/3 Ti 5/3 4 
is slightly larger than that for LiTi 2 C>4, as expected. A 
polynomial fit to the C P (T) of Li 4 / 3 Ti 5 / 3 04 above 12 K 
is shown by the dashed curves in Figs. § The uncer- 
tainties in the data and analyses for the Ti spinels have 
little effect on the analyses of C P (T) for LiV 2 04 in the 
following Sec. [II C , since as Figs. suggest, C lat (T) for 
LiV 2 04 is small compared to C e (T) of this compound at 
low temperatures. 

To quantify the difference above ~ 12 K between the 
C lat (T) of LiTi 2 4 and Li 4/3 Ti 5/3 04 noted above in 
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FIG. 2. Expanded plots below 20 K of the molar specific 
heat divided by temperature C p /T vs. temperature T of 
LiTi204 and IA4/3T1S/3O4, from Fig. |lj The solid curves 
are polynomical fits to the data for LiTi2 04, whereas the 
dotted curve is the inferred normal state behavior below 
T c — 11.8 K. The dashed curve is a polynomial fit to the 
data for LLj/aTis/aCU above 12 K. 



Sec. UlI A| , in Fig. |is plotted the difference AC lat (T) be- 
tween the measured C P (T) of Li 4/ 3Tig /3 04 and C lat (T) 
of LiTi 2 4 . The shape of AC lat (T) in Fig. | below 
~ 30 if is similar to that of an Einstein specific heat, 
but such a specific heat saturates to the Dulong-Petit 
limit at high T and does not decrease with T as the 
data do above 40 K. These observations suggest that 
intermediate-energy phonon modes in LiTi204 at some 
energy fc B T E2 split in Li 4 / 3 Ti 5/3 04 into higher (k B T E3 ) 
and lower (fcB?Ei) energy modes, resulting from the Li-Ti 
atomic disorder on the octahedral sites in Li4/ 3 Ti 5 / 3 04 
and/or from the difference in the metallic character of 
the two compounds. Following this interpretation, wc 
model the data in Fig. || as the difference ACg^stcin be- 
tween the Einstein heat capacities of two Einstein modes 
with Einstein temperatures of Tei and Tja. (neglecting 
the modes at high energy fc B TE 3 ), given byrJ 



AC]Jjnstein — 3rR 



Xl {T E1 /2T) 2 x 2 (T E2 /2T) 2 



sinh"(T E i/2T) sinh"(T E2 /2T) 



(4) 

where R is the molar gas constant, r — 7 atoms/formula 
unit and x\ and x 2 are the fractions of the total number 
of phonon modes shifted to Tei and away from Te 2 , re- 
spectively. A reasonable fit of the data by Eq. (ji|) was 
obtained with the parameters x\ = 0.012, Tei = 110K, 
X2 = 0.018 and T E 2 = 240 K; the fit is shown as the solid 
curve in Fig. |[ The model then predicts that a fraction 
(x 2 ~ x\)/x2 ~ 0.3 of the modes removed at energy &bTe2 
are moved to an energy /cbTe 3 3> /cbTe 2 . 
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FIG. 3. The difference AC lat between the lattice specific 
heats of Li4/3Ti 5 / 3 04 and LiTi20 4 vs. temperature T. The 
solid curve is a fit to the data by the difference between two 
Einstein specific heats in Eq. (Q) , whereas the dashed curve is 
the Schottky specific heat of a two-level system in Eq. (|B|). 
The error bars represent ±1% of the measured C P (T) for 

Ll4/3Ti5/304. 

An alternative parametrization of the experimental 
AC lat (T) data can be given in terms of the specific 
heat of |£u two-level system, described by the Schottky 
functioned 



^°Schottky 



xrR ( — 
\9i 



[1 + (So/5i)e 5 / T ] 



(5) 



where x is the atomic fraction of two-level sites, go and 
gi are respectively the degeneracies of the ground and 
excited levels and S is the energy level splitting in tem- 
perature units. Fitting Eq. (g) to the data in Fig. ||, 
we find gt/go = 4, x = 0.012 and 5 = 117 K. The fit is 
shown as the dashed curve in Fig.g. The accuracy of our 
AC lat (T) data is not sufficient to discriminate between 
the applicability of the Einstein and Schottky descrip- 
tions. 



5 



C. LiV 2 4 

Specific heat C P (T) data were obtained for samples 
from four batches of LiV2 04. Our first experiment was 
carried out on sample 2 (run 1) with mass 5 g. The C P (T) 
was found to be so large at low T (the first indication 
of heavy fermion behavior in this compound from these 
measurements) that the large thermal diffusivity limited 
our measurements to the 2.23-7.94 K temperature range. 
A smaller piece of sample 2 (0.48 g) was then measured 
(run 2) from 1.16 to 78. IK. Data for samples from two 
additional batches (sample 3 of mass 0.63 g, 1.17-29.3 K, 
and sample 4A of mass 0.49 g, 1. 16-39. 5K) were also 
obtained. Subsequent to the theoretical modeling of the 
data for sample 3 described below in Sec. [v], we obtained 
a complete data set from 1.14 to 108 K for sample 6 with 
mass 1.1 g from a fourth batch. A power series fit to the 
C P (T) data for sample 6 is shown as solid curves in Fig. [I]. 

We have seen above that C lat (T) of LiTi2C>4 is signifi- 
cantly different from that of Li4/ 3 Ti 5 / 3 04. Since L1V2O4 
is a metallic normal-spinel compound with cation oc- 
cupancies Li[V2]04 as in Li[Ti2]04, and since the for- 
mula weight of metallic LiTi204 is much closer to that 
of L1V2O4 than is that of the insulator Li4/ 3 Ti 5 / 3 04, we 
expect that the lattice dynamics and C lat (T) of L1V2O4 
are much better approximated by those of LiTi204 than 
of Li 4 / 3 Ti 5 / 3 04. Additionally, more precise and accu- 
rate C P (T) data were obtained for LiTi204 as compared 
to Li 4 / 3 Ti 5 / 3 04 because of the factor of three larger 
mass of the former compound measured than of the lat- 
ter. Therefore, we will assume in the following that the 
C lat (T) of LiV 2 4 from 0-108 K is identical with that 
given above for LiTi2C>4. We do not attempt to correct 
for the influence of the small formula weight difference 
of 3.5% between these two compounds on C lat (T); this 
difference would only be expected to shift the Debye tem- 
perature by < 1.8%, which is on the order of the accu- 
racy of the high temperature C P (T) data. The C e (T) of 
L1V2O4 is then obtained using Eq. (|2a|). 

The C e (T) data for samples 2 (run 2) and 6 of L1V2O4, 
obtained using Eqs. (||), are shown up to 108 K in plots of 
C* e (T) and C C {T)/T vs. T in Figs. |(a) and (b), respec- 
tively. An expanded plot of C C (T) below 9K for L1V2O4 
is shown in Fig. |](a), where data for sample 2 (run 1) 
and sample 3 are also included. The data for samples 2 
and 3 are seen to be in agreement to within about 1%. 
However, there is a small positive curvature in the data 
for sample 2 below ~ 3K, contrary to the small negative 
curvature for sample 3. This difference is interpreted to 
reflect the influence of the larger magnetic defect con- 
centration present in sample 2 as compared with that in 
sample 3, see Table B.El Therefore, we believe that the 
C e (T) data for sample 3 more closely reflect the intrinsic 
behavior of defect-free L1V2O4 compared to the data for 
sample 2 and all fits to C e (T) of L1V2O4 below 30 K by 
theoretical models to be presented in Sec. below are 
therefore done using the data for sample 3. As seen 
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FIG. 4. Electronic specific heat C e (a) and C e /T (b) vs. 
temperature T for LiV204 samples 2 (run 2) and 6. The 
error bars in (a) represent ±1% of the measured C P (T) for 
L1V2O4. 

in Fig. ||(a), the C C (T) data for sample 6 lie somewhat 
higher than the data for the other samples below about 
4 K but are comparable with those for the other samples 
at higher temperatures. This difference is also reflected 
in the magnetic susceptibilities x(T),S where x(T) for 
sample 6 is found to be slightly larger than those of other 
samples. 

To obtain extrapolations of the electronic specific heat 
to T = 0, the C C (T)/T data in Fig. | from 1 to 10 K for 
samples 3 and 6 were fitted by the polynomial 



^ = 7(0) + ^>„T 2 « , 



(6) 



71=1 



yielding 

7(0) = 426.7(6) mJ/molK 2 (sample 3) , (7a) 



7(0) = 438.3(5) mJ/molK 2 (sample 6) 



(7b) 
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FIG. 5. Expanded plot below 9K of the C c /T vs. T data 
for LiV2C>4 samples 2, 3 and 6. The solid and dashed curves 
are polynomial fits to the 1.1-10 K data for samples 3 and 6, 
respectively. 

The fits for samples 3 and 6 are respectively shown by 
solid and dashed curves in Fig. | The 7(0) values are an 
order of magnitude or more larger than typically obtained 
for transition metal compounds, and are about 23 times 
larger than found above for LiTi2C>4. 

The T-dependent electronic entropy S C (T) of LiV2 04 
was obtained by integrating the C e (T) /T data for sample 
6 in Fig. |(b) with T; the extrapolation of the C C {T)/T 
vs. T fit for sample 6 in Fig. [| from T = 1.16 K to 
T = yields an additional entropy of 5 e (1.16K) = 
0.505J/molK. The total S e (T) is shown up to 108K 
in Fig. ||; these data are nearly identical with those of 
sample 2 (run 2) up to the maximum measurement tem- 
perature of 78 K for that sample (not shown). The elec- 
tronic entropy at the higher temperatures is large. For 
example, if L1V2O4 were to be considered to be a strictly 
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sample 6 in Fig. fgb) with T. 
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FIG. 7. Measured specific heat divided by temperature 
Cp/T vs. temperature T for LiV204 sample 4A; correspond- 
ing data for sample 2 from Fig. [l] are shown for comparison. 
The lines are guides to the eye. 

localized moment system with one spin 5 = 1/2 per V 
atom, then the maximum electronic (spin) entropy would 
be 2i?ln(2), which is already reached at about 65 K as 
shown by comparison of the data with the horizontal 
dashed line in Fig. [| Our C P (T) data for one sample 
(sample 4A) of L1V2O4 were anomalous. These are shown 
in Fig. [?] along with those of sample 2 (run 2) for com- 
parison. Contrary to the C P (T)/T data for sample 2, the 
data for sample 4A show a strong upturn below ~ 5 K 
and a peak at about 29 K. We have previously associ- 
ated the first type of effect with significant (~ 1 mol%) 
concentrations of paramagnetic defects.la Indeed, Table | 
shows that this sample has by far the highest magnetic 
impurity concentration of all the samples we studied in 
detail. The anomalous peak at 29 K might be inferred 
to be due to small amounts of impurity phases (see Ta- 
ble ||). However, the excess entropy AS under the peak 
is rather large, AS" - 0.9 J/molK w 0.16i?ln(2). We also 
note that the height of the anomaly above "background" 
is at least an order of magnitude larger than would be an- 
ticipated due to a few percent of V4O7 or V5O9 impurity 
phases which order antiferromagnetically rwith Nccl tem- 
peratures of 33.3 and 28.8 K, respectively.cil It is possible 
that the 29 K anomaly is intrinsic to the spinel phase 
in this particular sample; in such a case Li-V antisite 
disorder and/or other types of crystalline defects would 
evidently be involved. As seen in Table Q, this sample 
has by far the largest room temperature lattice parame- 
ter of all the samples listed, which may be a reflection of 
a slightly different stoichiometry and/or defect distribu- 
tion or concentration from the other samples. Although 
these C P (T) data for sample 4A will not be discussed 
further in this paper, the origin of the anomaly at 29 K 
deserves further investigation. 
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FIG. 8. Linear thermal expansion coefficient a (left-hand 
scales) and a/T (right-hand scales) versus temperature T for 
L1V2O4 sample 6 from 4.4 to 297 K (a) and 4.4 to 50 K (b). 
The solid curves are the fit to the a(T) data by a polynomial. 



IV. THERMAL EXPANSION MEASUREMENTS 

The linear thermal expansion coefficient a(T) of 
L1V2O4 sample 6 was measured between 4.4 and 297 K. 
Figure ||(a) shows a(T) and a(T)/T over this T range, 
and Fig. ||(b) shows expanded plots below 50 K. At 
297K, a = 12.4 x 10~ 6 K" 1 , which may be com- 
pared with the value a w 15.6 x 10 _6 K _1 obtained 
for LiTi204 betweesu293 and 1073 K from x-ray diffrac- 
tion measurements.!^! Upon cooling from 297 K to about 
25 K, Qj-af L1V2O4 decreases as is typical of conventional 
metals.E3 However, a(T) nearly becomes negative with 
decreasing T at about 23 K. This trend is preempted 
upon further cooling below ~ 20 K, where both a(T) and 
a(T)/T exhibit strong increases. The strong increase i, 
a(T) below 20 K was first observed by Chmaissem et al. 
from high-resolution neutron diffraction data, which mo- 
tivated the present a(T) measurements. We fitted our 
ot(T) data by a polynomial in T over three contiguous 
temperature ranges and obtained the fit shown as the 
solid curves in Figs. ||. From the fit, we obtain 
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FIG. 9. Lattice parameter a a versus temperature T from 
4 to 297 K (a) and an expanded plot from 4 to 100 K (b) 
for LiV2 04. The filled circles are the neutroHndiffraction 
measurements of sample 5 by Chmaissem et al.Q'Q The solid 
curve is the linear thermal dilation obtained from our ca- 
pacitance dilatometer measurements of sample 6, assuming 
a o (0) = 8.22670 A. 

lim T ^o a{T) JT = 2.00 x 10~ 7 K~ 2 . 

Shown as the solid curve in Fig. ^|(a) is the linear ther- 
mal dilation expressed in terms of the lattice parame- 
ter a (T) = a o (0)[l + J a(T)dT], where we have used 
our polynomial fit to the a(T) data to compute a (T) 
and have set a o (0) = 8.22670 A. The a a (T) determined 
from thejieutron diffraction measurements by Chmais- 
sem et alE for a different sample (sample 5) are plotted as 
the filled circles in Fig. |[ The two data sets are in overall 
agreement, and both indicate a strong decrease in a (T) 
with decreasing T below 20 K. There are differences in 
detail between the two measurements at the lower tem- 
peratures as illustrated below 100 K in Fig. ||(b), suggest- 
ing a possible sample dependence. 

Our measurement of a(T)/T for sample 6 is com- 
pared with the measured C P (T)/T for the same sam- 
ple in Fig. |l^(a), where the temperature dependences 
of these two quantities are seen to be similar. We in- 
fer that the strong increase in a(T)/T with decreasing T 
below ~ 20 K is an electronic effect associated with the 
crossover to heavy fermion behavior. For most materials, 
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FIG. 10. (a) Comparison of the linear thermal expansion 
coefficient divided by temperature a(T)/T (left-hand-scale) 
with the measured specific heat C P (T)/T (right-hand-scale) 
for LiV204 sample 6. (b) Griineisen parameter F versus T, 
computed using Eq. ^ and the assumed value of the bulk 
modulus B given in the figure. 

the volume thermal expansivity (3 — 3a and C p are re- 
lated-through the dimensionless Griineisen parameter T, 
withEJ 



rc p 

B s Vm 



(8) 



where B s is the adiabatic bulk modulus and V^m is the 
molar volume. In this model, T = — dln^/dlnl^ where 
&(V) is a characteristic energy of the system. If indepen- 
dent contributions to C p can be identified, as assumed 
in Eq. (paj), a relation similar to Eq. ( paj ) exists for the 
thermal expansivity, with an independent T for each con- 
tribution 



lat/^lat 



p = & + /3 lat = 



T C C C + T M C, 
B s Vu 



(9) 



where C is understood to refer to measurements 
under constant pressure. For a metal, T c = 
d]nD*{E F )/d]nV = dln 7 (0)/dlnl/, and T lat = 
— d ln6*o / d InV. Here T>* (Ep) is the mass-enhanced quasi- 
particle density of states at the Fermi energy and the 
volume dependence of the elcctron-phonon interaction is 



neglected. Thus T c is a direct measure of the volume de- 
pendence of T>*(Ep). For a free electron gas, T c = 2/3. 
For most real metals T = ±3(2); e.g., T c = 0.92 ICu), 
1.6 (Au), 1.6 (V), -4.4 (Sr), -0.2 (Ba), 2.22 (Pd).IM 

We have computed T(T) for LiV 2 4 from Eq. ffl) us- 
ing the polynomial fit to our C p data for sample 6 and 
using the experimental a(T) data in Fig. || for this sam- 
ple. The molar volume of LiV 2 4 at low temperatures is 
given in Table [II. The bulk modulus is assumed to be 
B = 200(40) GPa, which is the range founco for the sim- 
ilar compounds Fe 2 3 , Fe 3 4 , FeTiQ*, MgO, Ti0 2 (rt 



tile), the spinel prototype MgAl 2 4 ,EJ and MgTi 2 5 .B 
The r obtained by substituting these values into Eq. (|^) 
is plotted versus temperature as the filled circles in 
Fig. |Io](b). Interpolation and extrapolation of T(T) is 
obtained from the polynomial fit to the a(T) data, shown 
by the solid curve in Fig. 0(b). From Fig. 0(b), T w 1.7 
at 108 K and decreases slowly with decreasing T, reaching 
a minimum of about 0.1 at 23 K. With further decrease 
in T, r shows a dramatic increase and wc obtain an ex- 
trapolated T(0) ps 11.4. A plot of T vs. T 2 obtained from 
our experimental data points is linear for T 2 < 30 K 2 , 
and extrapolates to 11.50 at T — 0, to be compared with 
11.45 as calculated from the smooth fitted relations for 
a(T) and C P (T); this justifies the (long) extrapolation 
of a(T) to T = 0. An accurate determination of the 
magnitude of T must await the results of bulk modulus 
measurements on LiV 2 04. Our estimated T(0) = r o (0) is 
intermediate between those of conventional nonmagnetic 
metals and those of /-electron heavy fermion compounds 
such as UPt 3 (r e = 71), UBei 3 (34) and CeCu 6 (571-with 
7(0) = 0.43, 0.78, and^67J/molK 2 , respectively.Eil 

From the expressionnJ relating C p to the specific heat 
at constant volume C v , and using our a(T) data and the 
estimate for B above, C V (T) of LiV 2 04 can be considered 
identical with our measured C P (T) to within both the 
precision and accuracy of our measurements up to 108 K. 



V. THEORETICAL MODELING: ELECTRONIC 
SPECIFIC HEAT OF LiV 2 4 

A. Single-Band Spin S — 1/2 Fermi Liquid 

As mentioned in Sec. [j| the high-temperature x(^) °f 
LiV 2 4 indicated a vanadium local moment with spin 
5 = 1/2 and g ~ 2. In the low-temperature Fermi liquid 
regime, for a Fermi liquid consisting of a single parabolic 
band of quasiparticles with S ri 1 /2 and iV G conduction 
electrons per unit volume 1/J23c3 the Fermi wavevector 
kp of LiV 2 4 assuming N c = 1.5 conduction electrons/ V 



atom is given in Table [II . In terms of the mass-enhanced 
density of states at the Fermi energy Ep for both spin 
directions r D*(Ep) 1 the 7(0) (neglecting electron-phonon 
interactions) and %(0) are given by 



7(0) 



] -V*(Ep) 



(10a) 



9 



x(o) 



,g 2 M | V*(E F ) 
4 1 + F n a 



(10b) 



where Fq is a Landau Fermi liquid parameter and 
1/(1 + Fq) = 1 — Aq is the Stoner enhancement fac- 
tor. The Fermi liquid scattering amplitudes A^ s are re- 
lated to the Landau parameters F^' s by Ag' s = F^' s /[l + 
F^' s /(2£ + 1)]. The superscripts "a" and "s" refer 
to spin-asymmetric an d sp in-symmetric interactions, re- 
spectively. Using Eq. ( 10a ) and the fcp value in Table Fjl, 
the experimental value of 7(0) for L1V2O4 in Eq. (7a[) 
yields the effective mass m* , Fermi velocity v-p , Fp , Fermi 
temperature Tp and T>*(Ep) for LiV2 04 giverwn Ta- 
ble |H From Eqs. (@) and ©, the Wilson ratioS R w is 
expressed as 



Rw 



1 



1 + F n a 



I -A* 



ID 



Substituting the experimental x(0-4-2K) = 0.0100(2) 
cm 3 /mol (Ref. § and 7(0) in Eq. @ for LiV 2 4 into 
Eq. (uh assuming g = 2 yields 



R w = 1.71(4) . 



(12) 



This i?w value is in the range of those found for 
many conventional as well as /-electron HF and IV 
compounds.?] The i?w value in Eq. (12) yields from 
Eq. ® 



F a 



-0.42, Al 



-0.71 



(13) 



In Fermi liquid theory, a temperature dependence is of- 
ten eamButed for C c at low temperatures having the 
forrrH e3 

TABLE III. Parameters for LiV20mri Abbreviations: for- 
mula weight FW; lattice parameter aopi3 (formula units/unit 
cell) Z\ theoretical mass density /9 calc ; molar volume 
Vm, itinerant electron concentration N c /V; Fermi wavevec- 
tor fcp = (37r 2 7V /V r ) 1//3 ; effective mass m* (free elec- 
tron mass m e ); Fermi velocity vf = hkF/m*; Fermi en- 
ergy Ef = h 2 fc|/2m* ; Fermi temperature Tp = Ff/&b; 
mass-enhanced density of states at Ff for both spin direc- 
tions, V"(E F ) = 3iV c /(2£;F) = m*fc P F/(7r 2 ft 2 ). 



Property 



value 



FW 

a (12K) 
Z 

p calc (12K) 
V M 

N c /V 
m* /m e 

I'F 

Ef 
T f 

V*(E f ) 



172.82 g/mol 
8.2269 A 

8 

4.123 g/cm 3 
41.92 cm 3 /mol 
4.310 xlO 22 cm -3 
1.0847 A" 1 
180.5 

6.96 x 10 5 cm/s 
24.83 meV 
288.1 K 
90.6 states/eV(V atom) 



C e (T) = T (0)T + ST 3 In (L ) + Oi T ' i 



(14) 



where 7(0) is given by Eq. (10a) and To is .-a. scaling 
or cutoff temperature. Engelbrecht and BedellEl consid- 
ered a model of a single-band Fermi liquid with the mi- 
croscopic constraint of a local (momentum-independent) 
self-energy, where the interactions are mediated by 
the quasiparticles themselves (in the small momentum- 
transfer limit). They find that only s-wave (I — 0) Fermi- 
liquid parameters can be nonzero and that the 5 coeffi- 
cient in Eq. ( O ) is 



3EB 



3tt 2 7(0) 



T 2 



24 



-A] 



< F a ' s < 00. 



(15) 



Within their 



where |^4g' s | < 1 and — | 
model, neither ferromagnetism nor phase-separation can 
occur. For Fq < 0, the only potential instability is to- 
wards antiferromagnetism and/or a metal-insulator tran- 
sition; in this case they find 1 < Fw < 2. For Fq > 0, a 
BCS superconducting state is possible and Fw < L The 
value of F a for LiV 2 4 in Eq. © is within the former 
range of this theory, ._. |—i|— 1 

Auerbach and Levirjfj and Millis et a^OQ formulated 
a Fermi-liquid theory of heavy electron compounds at low 
temperatures on the basis of a microscopic Kondo lattice 
model. The large enhancement of m* arises from the 
spin entropy of the electrons on the magnetic-ion sites 
(i.e., spin fluctuations).^ The Wilson ratio is Fw ~ 1.5 
and a T 3 lnT contribution to C e (T) is found. The ori- 
gin of this latter term is-Jiot ferromagnetic spin fluctu- 
ations ( "paramagnons" ) ,c3 but is rather electron den- 
sity fluctuations and the screened long-range Coulomb 
interaction.!! 3 ] The coefficient 5m of the T 3 In T term 
found by Milldfa is S M = n 2 k^V(l - 7r 2 /12)/5(^ F ) 3 , 
which may be rewritten as 



5M 



3tt 2 7(0) / 7T 



20T 2 



(16) 



Using the values 7(0) = 427mJ/molK 2 [Eq. ©], T F 
= 288 K (Table and A a in Eq. @, Eqs. @ and |TJ) 
respectively predict 



<5 EB = 0.0199 



mj 



molK 4 ' 



0.00135 



mj 
molK 4 



(17) 



We have fitted our low-temperature C e (T)/T data for 
LiV2C>4 sample 3 by the expression 

7 (T)^^=7(0) + <5T 2 lng)+eT 3 , (18) 

initially with e = 0. The fit parameters 7(0), 8 and To 
were found to depend on the fitting temperature range 
above 1 K chosen, and are sensitive to the precision of 
the data. The parameters obtained for 1-3 K and 1-5 K 
fits were nearly the same, but changed when the upper 
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FIG. 11. Electronic specific heat C c divided by tempera- 
ture T for LiV204 samples 2 (run 2) and 3 vs. T. The dashed 
curves are fits to the 1-5 K, 1-10 K and 1-15 K data for sam- 
ple 3 by the spin-fluctuation Fermi liquid model, Eq. ( |li| ) with 
e = 0, whereas the solid curve is a 1-30 K fit assuming £ / 0. 

limit to the fitting range was increased to 10 and 15 K. 
The fits for the 1-5 K, 1-10 K and 1-15 K fitting ranges 
are shown in Fig. along with the C C {T)/T data for 
sample 2 (run 2). As a check on the fitting parameters, 
we have also fitted the C C (T)/T data for sample 3 by 
Eq. (|l|) with s as an additional fitting parameter. The 
fit for the 1-30 K range is plotted as the solid curve in 
Fig. O. Since the fits for the smaller T ranges with e = 
and for the larger ranges with e ^ should give the most 
reliable parameters, we infer from the fit parameters for 
all ranges that the most likely values of the parameters 
and their error bars are 



7(0) = 428(2) 



mj 



mol K 2 ' 



1.9(3) 



mj 



mol K 4 



(19) 



The parameters in Eq. ( |l9| ) are very similar to those 
obtained using the same type of fit to C P (T)/T data 
for thei-heavy fermion superconductor UPt3 with T c = 
0.54 K,N for which^ylO) = 429-450 mJ/molK 2 and 
S = 1.99mJ/molK 4 .Ell!a Our T F and m*/m c values for 
LiV 2 4 (288 K and 181, Table |Tl|) are also respectively 
very similar to those of UPt 3 (289 K and 178).El 

The experimental 6 value in Eq. (|l^) is a factor of 
~ 10 2 larger than Seb and ~ 10 3 larger than 8m pre- 

-10 3 )] discrep-, 



dieted in Eq. fll7| ), A similar large [O(10 2 - 
ancy was found by Millis-for the 5 coefficient for UPts.cS 
As explained by MilliSjEj the large discrepancy between 
his theory and experiment may arise because the calcu- 
lations are for a single parabolic band, an assumption 
which may not be applicable to the real materials. How- 
ever, he viewed the most likely reason to be that his 
calculations omit some effect important to the thermo- 



dynamics such as antiferromagnetic spin fluctuations.0 
In this context, it is possible that the magnitude of 5 
predicted by one of the above two theories is correct, 
but that terms higher order in T not calculated by the 
theory are present which mask the T 3 hiX contribution 
over the temperature ranges of the fits;c3 in this case 
the large experimental 5 value would be an artifact of 
force- fitting the data by Eq. ( |l8[ ) . Indeed, we found that 
the fits were unstable, i.e., depended on the temperature 



range fitted (cf. Fig. 
of the theory of Milli 



). In addition, the applicability 
to LiV2 04 is cast into doubt by 
the prediction that the Knight shift at a nucleus of an 
atom within the conduction electron sea (not a "mag- 
netic" atom) "would be of the same order of magnitude 
as in a normal metal, and, would not show the mass en- 
hancement found in ^"cil In fact, the Knight shift of 
the 7 IajuiC|Lcus in L1V2O4 for T ~ 1.5-10K is about 
0.14%,HotaO which is about-6000 times larger than the 
magnitude (0.00024%) foundO at room temperature for 
the 7 Li Knight shift in LiTi 2 4 . Similarly, the 7 Li l/T x T 
from 1.5 to 4K in the highest-purity L1V2O4 samples is 
about 2.25 s -1 K -1 ,m3 which is about 600Q-,times larger 
than the value of 3.7 x lO^s^K" 1 foundO at 160K in 
LiTi204, where T\ is the 7 Li nuclear spin-lattice relax- 
ation time. 



B. Quantum-Disordered Antiferromagnetically 
Coupled Metal 

The antiferromagnetic (AF) Weiss temperature of 
LiV20 4 from x(T) measurements is \B\ — 30-60 K, yet 
the pure system exhibits neither static jaxitiferromagnetic 
AF nor spin-glass order above 0.02 K.qB A possible ex- 
planation is that the ground state is disordered due to 
quantum fluctuations. Wc consider here the predictions 
for C (T) of one such theory. A universal contribution to 
the temperature dependence of C e of a three-dimensional 
(3D) metal with a control parameter r near that required 
for a zero-temperature AF to quantum-disordered phase 
transition, corresponding to dynamical exponent z — 2, 
was calculated by Ziilicke and Millis,C3 which modifies the 
Fermi liquid prediction in Eq. (|l4|). Upon increasing T 
from T — in the quantum-disordered region, the system 
crosses over from the quantum disordered to a classical 
regime. The same scaling theory predicts that the low- 
T spin susceptibility is given by x(T) — x(0) + AT 3 /JL 
where the constant A is not determined by the theory.cSI 

Ziilicke and Millis found the electronic specific heat to 
be given by 24[ 



7o 



aRN, 



6T* 



F 



\rT* 



(20a) 



F(x) 



3V2 



dy^-2-Jl+yJl + 4x 2 y 2 . (20b) 
smh y 
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Here, 70 is the (nonuniversal) electronic specific heat co- 
efficent at T — in the usual Fermi liquid theory [7(0) 
above], T* is a characteristic temperature and Nq is the 
number of components of the bosonic order parameter 
which represents the ordering field: Nq = 3, 2, 1 for 
Heisenberg, XY and Ising symmetries, respectively. The 
number a is not determined by the scaling theory but is 
expected to be on the order of the number of conduction 
electrons per formula unit; thus for LiV2 04, we expect 
a ~ 3. We have defined F(x) such that F(0) = 1. The 
variable r is expected to be temperature dependent, but 
this temperature dependence cannot be evaluated with- 
out ascertaining the value of an additional parameter u 
in the theory from, e.g., measurements of the pressure 
dependence of C e (T); her e, we will assume r to be a 
constant .c3 From Eq. (20a), the T — value of 7 in the 
absence of quantum fluctuations is reduced by these fluc- 
tuations, and the measured 7(0) is 



7(0) = 7o 



aRN Qy 
6T* 



(21) 



We fitted our C c /T vs. T data for L1V2O4 sample 3 
by Eqs. (^0|), assuming Nq = 3. The fitting parameters 
were 70, a, r and T*; the 7(0) value is then obtained 
from Eq. (pl|). The 1-20 K and larger ranges did not 
give acceptable fits. The fits for the 1-5, 1-10 and 1- 
15 K fitting ranges are shown in Fig. [l2| From these fits, 
we infer the parameters and errors 



7« = 800(50)-^, 
T* = 18.9(4) K, 7(0) 



1 = 2.65(9) , r = 0.40(6) , 
430(1) mJ/mol K 2 . (22) 



Within the context of this theory, quantum fluctuations 
reduce the observed 7(0) by about a factor of two com- 
pared with the value 70 in the absence of these fluctua- 
tions. The value of a is close to the nominally expected 
value ~ 3 mentioned above. The relatively large value of 
r indicates that LIV2O4 is not very close to the quantum- 
critical point, and therefore predicts that long-range AF 
order will not be induced by small changes in external 
conditions (pressure) or composition. The former pre- 
diction cannot be checked yet because the required ex- 
periments under pressure have not yet been done. The 
latter expectation is consistent with the data available 
so far. Magnetic defect concentrations on the order of 
1% do induce static magnetic ordering below ~ 0.8 K, 
but this ordering is evidently of the short-range spin- 
glass typexJ Substitution of Zn for Li in Lii_. r Zn a ;V204 
induces spin-glass ordering for 0.2 < x < 0.9 but long* 
range AF ordering does not occur until 0.9 < x < 1.0.0 
Finally, two caveats regarding the fits and discussion in 
this section are in order. The first is that (unknown) cor- 
rections of-order (T/T*) 2 and r 1 to the theory of Ziilicke 
and MillisHa exist but have not been included in the pre- 
diction in Eqs. (pOj); incorporating these corrections may 
alter rtiie parameters obtained from fits to experimental 
data£3 The second caveat is that the theory 
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FIG. 12. Electronic specific heat divided by temperature 
C /T vs. T for L1V2O4 samples 2 (run 2) and 3. Fits 
the data for sample 3 by the theory of Ziilicke and Millis,Ej 
Eqs. ([20]), are shown for the fitting ranges 1-5 K (long-dashed 
curve), 1-10 K (solid curve) and 1-15 K (short-dashed curve). 

may need modification for compounds such as LiV2 04 in 
which geometric frustration for AF ordering exists in the 
structured 



C. Spin-1/2 Kondo Model 

Calculations of the impurity spin susceptibility x(T) 
and/or impurity electronic contribution C e (T) to the spe- 
cific heat for the S — 112. Kondo model were carried out 
by Wilsonfl and othersE^TEa using different techniques. 
Both x(T) and C C (T) depend only on the scaling param- 
eter T/Tk, where Tk is the Kondo temperature (here, 
we use Wilson's definition! 2 ]) . The impurity x(T) is pre- 
dicted to be Curie- Weiss-like at temperatures high com- 
pared with Tk , and to level out at a constant high value 
for T < 0.1 Tk due to the formation of a singlet ground 
state. 1 — - 

In the limit of zero temperature, one hasE3 



7 (T - 0) 



irWNk B 
6T K 



(23) 



where M is the numbep-n£,impurity spins. The Wilson 
numbeia W is given byE3 ! E2l 



W 



1/4 - 
7e ' 7r 



1/2 



1.290 268 998 



(24) 



where In 7 w 0.577 215 664 902 is Euler's constant. Set- 
ting N = Na, Avogadro's number, one obtains from 
Eqs. ( p3| ) and ( pi| ) the electronic specific heat coefficient 
per mole of impurities 
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7(0) = 



ttWR _ 5.61714 J/molK 



(25) 



Fits: S = 1/2 Kondo Model 



6T K T k 

To characterize the T dependence of C e , we utilized 
accurate numerical calculations using the Bethe ansatz 
by Jerez and Andrei. E3 The calculated C e (T) shows a 
maximum, m&x[C c (T)/Nk-B] — 0.177275, which occurs 
at T max /T K = 0.6928. The calculations were fitted by 
the expressions 

Ce(T) 



Nkv 



/(*) 



9(t) = — 



/(*) 



C e {T) 
Nk B T/T K 

t(l + ait + a 2 t 2 + a 3 t 3 + a 4 i 4 ) 



(26a) 



(26b) 



irW 
~6~ 



1 + a 5 t + aat 2 + a 7 t 3 + a s t 4 + a 9 t 5 



(26c) 

where t = T/Tk and the coefficients a n for the two 
types of fits are given in Table IV for the fitting range 
0.001 < t < 100. Equations (|26j) incorporate the zero- 
temperature limit in Eqs. (123 |25| ). The maximum (rms) 
deviations of the C e (T) fit from the calculated numerical 
data are 0.011% (0.0035%) for < t < 3 and 0.031% 
(0.021%) for 3 < t < 10 but then progressively deterio- 
rate to 0.48% (0.14%) in the region 10 < t < 92. The 
corresponding deviations for the C e (T)/T fit are 0.0044% 
(0.00091%), 0.031% (0.017%) and 5.1% (1.6%). 

The experimental C C (T)/T data for L1V2O4 sam ple 3 
were least^aquare fitted from 1.2 to 5K by Eqs. (|26l 
and (p6c|)j£j yielding T K , and then 7(0) from Eq. (|2£ 




T K = 26.4(1) K , 7 (0) = 426(2) mJ/molK 2 . 

The fit is shown in Fig. [l3] as the solid curves. For com- 
parison, also shown in Fig. |l^(a) are the predictions for 
Tk = 25 K and 28 K. Unfortunately, despite the good 
agreement of the theory for Tk = 26.4 K with our mea- 
sured C C (T) at low T, the S — 1/2 Kondo model predic- 
tion for x(T) qualitatively disagrees with the observed 
temperature dependence at low TB This difficulty of self- 
consistcntly fitting the C C (T) and x(T) data is a problem 



TABLE IV. Coefficients a n in Eq. (|26c|) in the fits to the 
theoretical prediction for the specific heat vs. temperature of 
the 5 = 1/2 Kondo model by Jerez and Andrei .E3 





C(T) Fit 


C{T)/T Fit 


ai 


9.1103933 


6.8135534 


(12 


30.541094 


21.718636 


az 


2.1041608 


2.3491812 


0,4 


0.0090613513 


0.017533911 


as 


9.1164094 


6.8158433 


(l 6 


36.143206 


27.663307 


at 


67.91795 


48.229552 


as 


53.509135 


40.216156 


n» 


1.7964377 


2.4863342 



o 
E 
—> 
E 



O 




150 
100, 



LiVX> 

2 4 

Sample #3 
Sample #2 (Run #2) 
T K = 26.4 K 

T K = 25.0 K 

T K = 28.0 K 



a i 



10 15 20 25 30 
T(K) 




10 20 30 40 50 60 70 80 
T(K) 

FIG. 13. (a) Electronic specific heat divided by tempera- 
ture C e (T)/T data for LfV^CU samples 2 (run 2) and 3 below 
30 K (open symbols) and a fit (solid curve) of t he dat a fo r 
sample 3 by the 3 = 1/2 Kondo model, Eqs. ( p6b| ) and ( p6c| ), 
for a Kondo temperature Tk = 26.4 K. Shown for comparison 
are the predictions for Tk = 25.0 K (long-dashed curve) and 
28.0 K (short-dashed curve), (b) The same data and the fit 
with T K = 26.4 K in a plot of C c vs. T up to 80 K. 

we have encountered in all our attempts so far to fit our 
data for both measurements over any extended tempera- 
ture range by existing theory (see also the next section). 



D. Local Moment High- Temperature Description 

As discussed above, the x(T) data for L1V2O4 sug- 
gest that at high temperatures a local moment descrip- 
tion in which the moments are antiferromagnetically cou- 
pled with Wfeiss temperature 9 ~ —30 to — 60 K may be 
applicable. Accordingly, we have calculated the mag- 
netic specific heat C m (T) for localized moments on the 
octahedral (B) sublattice of the ^4[i?2]04 spinel struc- 
ture assuming nearest-neighbor AF Heisenberg interac- 
tions using the general high-temperature series expansion 
(HTSE) results of Rushbrooke and Wood.E3 The Hamil- 
tonian is H — J J2<ij> &i ' ^31 where the sum is over all 
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exchange bonds and the exchange constant J > corre- 
sponds to AF interactions. In terms of this Hamiltonian, 
= —zJS(S + l)/3, where z = 6 is the coordination 
number for the B sublattice of the spinel structure. The 
above range of 9 then gives J/k-g = 20-ADK assuming 
5 = 1/2. The general HTSE prediction isEI 



C m (T) = 



6t 2 



E 



n(S) 
t n 



(28) 



where t = k^T / J and the coefficients c n depend in gen- 
eral on the spin-lattice structure in addition to S. The 
coefficients c n for the B sublattice of the spinel structure 
with S = 1/2 and S = 1 up to the maximum available 
n max _ g are gj ven m Table fv[ The predictions for C m 
versus scaled-temperature kBT/[JS(S + 1)] with n max = 
5 are very similar for S — 1/2 and S = 1. A comparison 



of the C m (T) predictions for n max — to 5 indicates 
that the calculations for n max = 5 are accurate for 
k B T/[JS(S + 1)] > 2.5, a T range with a lower limit 
slightly above the temperatures at which broad maxima 
occur. 

In Fig. |lj the HTSE prediction of C m (T) for the B 
sublattice of the spinel structure with n max = 5, S = 1/2 
and J/fcs = 20 K in Eq. ( |28| ) is compared with the ex- 
perimental C C (T) data for LiV204 samples 2 and 6 from 
Fig. |(a). The HTSE C m (T) has a much lower magni- 
tude than the data and a qualitatively different temper- 
ature dependence. From Eq. ( p8[ ), changing J just scales 
the curve with T . Thus the local moment picture is in 
severe disagreement with our C e (T) measurements, de- 
spite the excellent agreement between the corresponding 
HTSE xffA prediction and the X {T) data from 50-100 K 
to 400 KM 
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FIG. 14. Comparison of the high temperature series expan- 
sion prediction for the magnetic specific heat C m (T) of the B 
sublattice of the A[i?2]04 spinel structure assuming S = 1/2, 
J/k B = 20 K and n max = 5, given by Eq. @ with c n co- 
efficients in Table [v], with the experimental C C (T) data for 
LiV2 04 sample 2 (run 2) and sample 6 from Fig. KM a). 

TABLE V. Coefficients c„ in Eq. @ for the high temper- 
ature series expansion for the magnetic specific heat of the B 
sublattice of the spinel structure, for the indicated values of 
spin S. 



n S= 1/2 S = 1 

1 -1/2 -13/6 

2 -23/16 -3 

3 65/48 715/36 

4 1183/768 -4421/324 

5 -18971/7680 -670741/6480 



VI. SUMMARY AND CONCLUDING REMARKS 

We have presented C P (T) data for LiVp04 sample 6 
which extend our previous measurementsa up to 108 K. 
We have also presented C P (T) data for the isostructural 
superconducting compound LiTi204 (T c = 1L8K) up 
to 108 K which complement our earlier datata on the 
isostructural nonmagnetic insulator Li4/ 3 Ti 5 / 3 04. We 
concluded here that the lattice contribution C lat (T) to 
C P (T) for LiTi 2 04 provides the most reliable estimate 
of the C lat (T) for LiV2C>4, and we then extracted the 
electronic contribution C C (T) to C P (T) of LiV204 from 
1.2 to 108 K. Inelastic neutron scattering measurements 
of the lattice dynamics and spin excitations would be 
very useful in interpreting the measurements presented 
here. It will be important to determine whether or not 
there exist significant differences in the lattice dynamics 
of LiV2C>4 and LiTi 2 04; in our data analyses and model- 
ing, we have assumed that these compounds are similar 
in this respect. 

For two high-magnetic-purity LiV2 04 samples 3 and 6, 
the electronic specific heat coefficients 7(T) = C C (T)/T 
were found to be 7(1 K) = 0.42 and 0.43 J/moleK 2 , re- 
spectively. To our knowledge, these values are signif- 
icantly larger than previouslyp&eported for any metal- 
lic transition metal compoundEJ For LiTi204, we found 
7 = 0.018 J/moleK 2 . 7 (T) of LiV 2 4 decreases rapidly 
with increasing temperature from 4 to 30 K and then de- 
creases much more slowly from a value of 0.13 J/molK 2 
at 30 K to 0.08 J/molK 2 at 108 K. Even these latter two 
7 values are exceptionally large for a metallic d-electron 
compound. The temperature dependences of 7, x> t ne 
low-T resistivity and the 7 Li NMR properties are re- 
markably similar to those of the heaviest mass /-electron 
heavy fermion compounds!!! In a plot of x(0) versus 7(0), 
the data point for LiV2C>4 sits amid the cluster formed 
by the /-electron heavy fermionparid intermediate valent 
compounds as shown in Fig. |15|,lj where several 
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FIG. 15. Log-log plot of the magnetic susceptibility x(0) 
versus electronic specific heat coefficient 7(0) at zero tem- 
perature for a variety of /-electron heavy fermion and in- 
termediate valence compounds compiled from the literature 
(after Ref.|^). The plot also includes data for several elemen- 
tal and/or d-electron metals and our data point for LiV2 04. 
Here, a "mol" in the axis labels refers to a mole of transition 
metal atoms for the d-metal compounds, and to a mole of 
/-electron atoms for compounds containing lanthanide or ac- 
tinide atoms. The straight line corresponds to a Wilson ratio 
Rw = 1 for quasiparticles with spin 5 = 1/2 and g- factor 
g — 2, which is also the Wilson ratio for a free-electron Fermi 
gas. 

data for elemental metals, the A-15 superconductor VsSi 
(T c = 17K),L30 and superconducting ^and/or metallic 
d-metaL-oxides LiTi 2 4 (T c < I3.7K),EJ Sr 2 Ru0 4 (T c 
= 1K)J23 and (Vo.95Tio.o5)203,Ej are also included for 
comparison. 

From our theoretical modeling in Sec. [v|, Fermi liquid 
models and the S — 1/2 Kondo model (with a Fermi 
liquid ground state) are capable of describing our C e (T) 
data for L1V2O4 from 1 K up to ~ 10 K, although the 
magnitudes of the derived parameters remain to be un- 
derstood theoretically. The localized moment model in 
Sec. VD failed both qualitatively and quantitatively to 
describe the data. None of the models we used can ac- 
count for the additional contribution to C C (T) at higher 
temperatures, from ~ 10 K up to our high temperature 
limit of 108 K, which appears to be distinct from the con- 
tributioa,beginning at much lowe r T a nd could arise from 
orbital,EaE!J charge and/or spinOO excitations. The 
crystalline electric field and/or the spin-orbit interaction 
may produce some energy level structure whichds ther- 
mally accessible within our temperature range.L3 Con- 
ventional band structure effects cannot give rise to our 



results. 

As is well-known for conventional metals, the electron- 
phonon interaction increases 7 by the factor (1+A), where 
A is the electron-phonon coupling constant, but do es no t 
affect x, i-e-, V*(E F ) -> V*(E F )(1 + A) in Eq. (jlO^). 
One can correct the observed Wilson ratio for electron- 
phonon interactions by multiplying the observed value by 
(1 + \)L3 The electron-phonon interaction is not taken 
into account in any of the analyses or modeling we have 
done. This correction would have had a significant quan- 
titative impact on our analyses if we used, e.g. A ~ 0.7 
as in LiTi204 (Rcfs. pfl| , p8| ); most previous analyses of 
the specific heats of other (/-electron) HF compounds 
also did_not take the electron-phonon interaction into 
account B 

From our combined specific heat and thermal expan- 
sion measurements on the same sample 6 of LiV204 from 
4.4 to 108 K, we derived the Griineisen parameter T(T) 
which shows a dramatic enhancement below ~ 25 K as 
the compound crosses over from the quasilocal moment 
behavior at high temperatures to the low-temperature 
Fermi liquid regime, confirming the discovery of Chmais- 
sem et al. from neutron diffraction measurements. B Our 
estimated extrapolated value of the electronic Griineisen 
parameter r e (0) is about 11.4, which is intermediate be- 
tween values for conventional metals and for /-electron 
heavy fermion compounds. This large value indicates 
a much stronger dependence of the mass-enhanced den- 
sity of states on the volume of the system than simply 
due to the decrease in the Fermi energy with increasing 
volume as in the quasi-free electron picture. In the /- 
electron HF systems, the large r o (0) values are thought 
to arise from a strongly volume dependent hybridization 
of the /-electron orbitals with those of the conduction 
electrons ESO In the present case of LiV20 4 , the origin 
of the large r e (0) is unclear. 

It is conceivable that the same mechanism is respon- 
sible for the heavy fermion behavior in LiV204 as in 
the /-electron heavy fermion systems if one of the 1.5 d- 
electrons/V atom is localized on each V atom due to elec- 
tron correlation effects and-crystalline electric field or- 
bital energy level structure,!!] and if the orbital occupied 
by the localized electron is hybridized only weakly with 
the conduction electron states. That such localization 
can occur in similar systems is-s.upported by calculations 
for the d 1 compound NaTi02-II3 Additional scenarios for 
the heavy fermion behavior mechanism(s) are given by 
Kondo et alu U involving the geometric frustration for AF 
ordering within the V sublattice and/or low- lying cou- 
pled dynamical orbital-charge-spin excitations. Further 
experimental and theoretical investigations of the physi- 
cal properties of LiV204 may thus reveal interesting new 
physics which may also allow a deeper understanding of 
the /-electron heavy fermion class of materials. 
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